Determination of 817(0) from published data 
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The experimental landscape for the ^Be+p radiative capture reaction is rapidly changing as new 
high precision data become available. We present an evaluation of existing data, detailing the 
treatment of systematic errors and discrepancies, and show how they constrain the astrophysical 
S factor (Sit), independent of any nuclear structure model. With theoretical models robustly 
determining the behavior of the sub-threshold pole, the extrapolation error can be reduced and a 
constraint placed on the slope of Sn. Using only radiative capture data, we find 517(0) = 20.7 ± 
0.6(stat) ± l.Q[syst) eV b if data sets are completely independent, while if data sets are completely 
correlated we find Si7(0) = 21.4 ± 0.5(siat) ± lA{syst) eV b. The truth likely lies somewhere in 
between these two limits. Although we employ a formalism capable of treating discrepant data, 
we note that the central value of the S factor is dominated by the recent high precision data of 
Junghans ei al., which imply a substantially higher value than other radiative capture and indirect 
measurements. Therefore we conclude that further progress will require new high precision data 
with a detailed error budget. 

PACS numbers: 25.70.De, 26.20.+f, 26.65. +t, 27.20. +n, 07.05.Kf 



I. INTRODUCTION 

The comparison of measured and predicted ®B solar 
neutrino fluxes represents a test of solar models and an 
opportunity to learn about the properties of the elec- 
tron neutrinos produced in the Sun. Recent measure- 
ments at SNO 1] have determined the total flux of ac- 
tive neutrinos emitted in the Z?"*" decay of ^B with a com- 
bined statistical and systematic precision of 9%. The 
implications of solar and reactor neutrino flux measure- 
ments for neutrino mixing parameters are explored in, 
e-g-, Q;y|' Theoretical predictions of the *B solar neu- 
trino flux are now substantially more uncertain than the 
experimental measurements. Recent predictions have un- 
certainties of 20% j4| and 15% Q. The latest effort 
to estimate the theoretical uncertainty found a value of 
23% 0. This error is completely dominated by the un- 
certainty in the heavy element abundance of the Sun, 
which has recently been revised to a level 2.5 times larger 
than the previous adopted value Q- The contribution 
to the theoretical error budget made by S17 has now 
been estimated at 3.6% based on the recommendation 
of Si7(0) = 21.4 ± 0.5 (expt.) ± 0.6 (theor.) eV b given 
in Q. Even if the uncertainty due to S17 represents a 
small fraction of the total theoretical uncertainty, in our 
view an independent determination of its value based on 
the available experimental data and theoretical models is 
worthwhile. Here we provide reliable determinations of 
the total uncertainty in the S factor, and flnd that nearly 
all previous analyses have underestimated the error. 

In this paper we detail our fitting procedure, presenting 
formalisms for propagating systematic errors and com- 
bining multiple data sets, some of which were first pre- 
sented in Q. We then describe the available data rele- 
vant for ^Be(p,7)®B, and present the data sets used in 
this analysis. Next we briefly discuss the structure mod- 
els used to extrapolate experimental data to solar ener- 



gies, and the question of how the models can be tested. 
Finally, we present our constraints on the astrophysical 
S factor using: (1) the energy dependence of the best 
structure models, (2) a pole model parametrization inde- 
pendent of structure models, and (3) a constrained pole 
model parametrization, where theory is used to robustly 
determine the behavior of the subthreshold pole. 



II. FITTING PROCEDURE 

It is often desirable to determine an average or best-flt 
representation of experimental data, and the uncertain- 
ties in such a representation. The standard techniques 
are generally not discussed in the literature, and are as- 
sumed to be well known. We detail here a formalism 
needed to properly take into account correlated data and 
systematic errors. 

When modeling data, one generally uses a maximum 
likelihood or minimum y^ formalism to determine the 
best fit. The x^ is usually defined as: 



E 



V{xi) - Hi 



(1) 



where fii, and Ui are the mean and standard deviation of 
the i*'* data point and y{x) is the value calculated from 
the model we are using to describe the data as a function 
of the independent variable x (e.g. y{x) -^ Sn{E)). As 
an example, a linear model consists of P parameters and 
P basis functions (i.e. y{x) = '^„^iO,pYp{x)). Minimiz- 
ing x^ defined in this way also minimizes the variance in 
the best fit model, but is limited in that it assumes the 
data points are independent and there is no explicit pre- 
scription for dealing with systematic errors or multiple 
data sets. A more general treatment has been described 
recently in an analysis of reaction cross sections relevant 



for big bang nucleosynthesis y\. We shall adopt a sim- 
ilar procedure for studying the ^Be(p,7)*B reaction. A 
clear understanding of how statistical and systematic un- 
certainties propagate through the analysis is an essential 
aspect of this approach. It is assumed that the dominant 
systematic error is the normalization uncertainty, which 
can be parametrized by a relative uncertainty, e„, for 
each data set n. This global normalization error induces 
correlations between data points in the same data set. 
The correlation between two data points is given by [3 : 



C,. 



= (1 + <^l)<^fJ^„J,^ + 4M^„AiJ„ 



(2) 



where ct^^ is now the statistical uncertainty of the i*^ data 
point of data set n and 5i„,j„ is the Kronecker delta. Us- 
ing this correlation matrix, one can define a more general 
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where we have left off the data set subscript n for clarity. 
Given the simple form of the covariance matrix in Eq. |21 
the inverse can be found analytically |9j, and we obtain 
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where we have again suppressed the index n. When the 
systematic errors are smaller than the statistical errors 
(e„/ii^ < ai^) or in the limit of large data sets, the x^ of 
Eq. 131 reduces to that of Eq. ^ with the statistical error 
weighting the x^ rather than the total error. 

We stress that the x^ quantity is a statistical device 
only; its minimum value tells us only how good the fit is 
within statistical uncertainties and its curvature tells us 
only about the statistical uncertainties in the fit. This is 
seen as a 1/VN scaling in the best fit uncertainty, where 
N is the number of data points. Systematic errors are 
not reduced by this factor. The total uncertainty then 
consists of the statistical error and the intrinsic normal- 
ization error. Statistical errors contain no information 
about the quality of the fit. This can be seen by arbi- 
trarily shifting data points away from each other, keeping 
their absolute uncertainties the same. We need an addi- 
tional measure of how far data fall from the best fit. How 
then do we address the quality of the fit? 

In Ref. [3, such a discrepancy error is defined as a 
measure of the fit quality. The discrepancy error is the 
weighted dispersion of the data relative to the best fit: 



T,t,jC,j[yixt) ~ ^i^][y{xj) - ptj] 
J2r.TC~jy{xt)yixj) 



(5) 



The absolute size of this discrepancy error tells us how 
well the data are described by the best fit, while its 
size relative to the intrinsic normalization error quantifies 
possible unknown systematics. One may be tempted to 



reduce the size of the discrepancy error by the number of 
degrees of freedom (e.g. x^ ~^ ^ 1^)^ but this is inappro- 
priate because it assumes that the unknown errors we are 
trying to take into account can be propagated through 
the data analysis. As discussed earlier, systematics are 
not reduced by \/\/N as are statistical errors. We adopt 
a total normalization error defined as the quadrature sum 
of the intrinsic normalization error and this discrepancy 
error, as was done in Ref. :9j|. 

We summarize our procedure for analyzing single data 
sets as follows. 

1. We find best fits and statistical uncertainties, where 
the statistical errors of the data points dominate 
the x^ analysis. 

2. The total normalization error is the quadratic sum 
of the intrinsic normalization error and our quality- 
of-fit measure, the discrepancy error. 

We are then left with the remaining question of how to 
treat multiple data sets. We discuss two methods. The 
first method, adopted in p treats the data sets as totally 
correlated and comprising a single data set. The second 
method treats the data sets as if they were completely 
independent. Reality is likely between these two possi- 
bilities; however at present we have no good prescription 
for determining how correlated data sets actually are, and 
therefore present results for these two limiting cases. Ad- 
ditionally, both methods yield similar results, suggesting 
that the methods are accurate and robust. Generally, the 
totally correlated method yields more conservative un- 
certainties as the intrinsic normalizations arc not treated 
statistically. 



A. Completely Correlated Data Sets 

Assuming that individual data sets are totally corre- 
lated simplifies the analysis. Since the exact nature of 
the correlation is unknown, we cannot rigorously define 
a correlation matrix, so we continue to use the definition 



jm — Sn,mCi„,jn) and its inverse in Eq. 0] 
%C~^, ), generalized for multiple data sets. 



in Eq. (a,„ 

By virtue of the large data set limit for the inverse covari- 
ance matrix, we are led to the conclusion that the precise 
nature of the correlations is relatively unimportant, as 
the statistical uncertainties dominate the x^ analysis. 

As before, the total normalization error is the 
quadratic sum of the intrinsic normalization error and 
the discrepancy error (e^^^ = e^^^^ -^ e^^^j. The dis- 
crepancy error as defined in Eq. [S] is valid for a single 
data set and hence also for the case of totally correlated 
data sets. Since the individual data sets are completely 
correlated, the overall normalization error must be some 
average of the individual normalization errors. We adopt 



the normalization error prescription of Ref. |9j : 



(6) 



where e„ are the individual data set normalization er- 
rors, and Xn is the x^ P^r datum of data set n with 
respect to the best fit. This weighting scheme gives more 
weight to data sets that agree with the best fit model. We 
also point out that this normalization error assignment 
is bounded by the smallest and largest normalization er- 
rors and is not reduced by the number of data sets. This 
reflects the fact that the data sets are completely cor- 
related and the normalization errors cannot be treated 
statistically in this case. 

The expectation value and statistical variance of the 
best fit model are denoted E[y{x)] and V5[2/(x)], re- 
spectively. Including the normalization error, the to- 
tal variance in the best fit model is VT[2/(a;)] = (1 + 
e'^gf)Vs[y{x)] + €fgfE[y{x)]'^ . This completes our descrip- 
tion of the formalism for completely correlated data sets, 
which we refer to as the correlated normalization analy- 
sis. 



B. Completely Independent Data Sets 

If data sets are truly independent from each other, we 
can treat the normalizations statistically, once best fits 
and statistical errors are found for each data set. We con- 
sider only linear models, for which y{x) = X]p=i flpYp(a;), 
where the Yp{x) are known functions, and follow the 
prescription laid out in the previous section for each 
data set. The expectation value and statistical variance 
of y{x) are E[y{x)] = Y.p=i^P^p{^) and Vs[y{x)\ = 
Spg=i^p,gYp(a:)Yq(a;), respectively. With the best fit 

parameters CLp and parameter covariance C^q for each 
data set n in hand, we can combine individual data sets 
and find a global best fit. Upon minimization, the x^ can 
be decomposed as: 



x' = xL„ + E E c^rr\«p - 4"))(«, - 4")), (7) 

n p,g— 1 

where the calligraphic C^q is the covariance between the 
p*^ and q*"^ parameters of the n*-^ data set, with best fit 
parameters given by Op , provided there are at least as 
many data points N as fitting parameters P. Note that 
with linear models all parameters are gaussian, i.e., the 
X^ is a quadratic function of the parameters. 

In order to combine multiple data sets, we must first 
replace the statistical variance, Vs[y(a:)], with the total 
variance VT[2/(a;)]. In other words, we replace the pa- 

(n) 

rameter covariance matrix Cp^q for each data set n with 
(l-|-e(o( „)Cp"q -|-e(oj_„ap Sg . Inserting this new param- 
eter covariance into Eq. [7| we minimize the x^ to find 



the best fit parameters and their statistical covariances, 
which now include the normalization errors of each data 
set. 

We also need to quantify how well data sets agree with 
each other. The variance of the best global fit contains no 
information about the mutual consistency of data sets. 
This is seen if we shift the mean values of data points 
in a single data set, but keep the same absolute uncer- 
tainties. The parameter covariance is not changed, even 
though data sets can be severely discrepant. Thus we 
need to define a global discrepancy error. To do this, we 
calculate the renormalizations needed for the global best 
fit to minimize the x^ for individual data sets. The dis- 
persion of these renormalizations provides an estimate of 
the discrepancy between data sets. We find renormaliza- 
tions a„ and their total errors (t„ for each data set and 
then their dispersion according to 



E„ 



2-^ix 



(8) 



This definition is similar to that discussed for the in- 
dividual data sets, except that here no correlations ex- 
ist between data sets. Given the global best fit E\iy{x)\ 
and its statistical variance Vs\y(x)\ the total variance in 
this case is given by: VT\y{x)\ = (1 -h eLc)^s[y(a;)] + 
^disc^\n^^)^ ■ We refer to this formalism as the indepen- 
dent normalization analysis. 



III. THE DATA 

With these formalisms in place, we now discuss the 
data available to constrain the astrophysical S factor of 
the ^Be(p,7)*B reaction. In some cases there is suffi- 
cient reason to exclude data sets from the analysis. We 
consider only low energy data, Ecm < 425 keV, when 
determining the best fits to Sij{E) as nuclear struc- 
ture uncertainties complicate and render more uncertain 
the ex trap olation when higher energy data are included 

sua III 



A. Radiative Capture Data 

Initialh^ we considered the data sets of Kavanagh 12 1, 
Parker [13, Vaughn et al. [lj|, Filippone et al. ^|, 
Strieder et al. JI6l| . Hammache- 1 Il7| . Hammache-2 '18], 
Hass et al. [I^, Junghans-BEI |23|, Junghans-BE3 ^, 
and Baby et al. [211. We exclude the data of Ka- 
vanagh [131, Parker uM, and Vaughn et al. [lj| because 
these authors do not present enough information to ade- 
quately determine a normalization error. We do not use 
the measurement of Hass et al. [I^ simply because the 
data lie above our 425 keV energy cutoff. 

Several details of our analysis bear mention here. One 
of the two target thickness determinations in the Filip- 



pone et al. measurement [T^ relies on the ^Li(d,p) reac- 
tion. We adopt the recommendation of |23| for the value 
of this reaction cross section. The Hammache-2 [l^ data 
consist of 3 points, two of which are measured relative to 
the third. Ideally, one would like to include all 3 points, 
but not enough information is given on the third point 
to determine an intrinsic normalization error. We thus 
adopt the 2 relative measurements as the data set, using 
the third to determine the normalization error. Ref. [g 
presents data from their BEl measurement renormalized 
using the BE3 data. These renormalized BEl data are 
not independent of the BE3 data. Therefore, we consider 
here the BE3 data ^ and the original BEl data |20| . 
which are independent. Finally we note that there is 
some discussion in the literature |3] r egar ding the uncer- 
tainties in the data of Baby et al. W\\. We take the 
uncertainties for this measurement from Table II of Ref. 

m. 



B. Coulomb Dissociation Data 

We considered the Coulomb dissociation (CD) data of 
Kikuchi et al. [23, |2j, Iwasa et al. [23, Schiimann et 
al. [2g and Davids et al. IllL l27l . We exclude the mea- 
surements of Kikuchi et al. [23, |2J| and Iwasa et al. [23 
due to concerns over the way these data were analyzed. 
The data from these measurements were not analyzed 
using a cut on the maximum scattering angle of the ^B 
center-of-mass, corresponding semiclassically to a mini- 
mum impact parameter. This means that the effects of 
nuclear absorption, diffraction dissociation, and E2 tran- 
sitions are present in the data, and that the inferred E\ S 
factors may not be reliable. Since we lack the detailed 
experimental information required to correct for these 
effects, we do not consider these data here. Hence we 
include only the CD data of Schiimann et al. |2g and 
Davids et al. [I^, 1231 , which were analyzed using scatter- 
ing angle cuts to minimize the nuclear and E2 contribu- 
tions and their associated uncertainties. 



IV. THEORETICAL STRUCTURE MODELS 

In general, two classes of models have been used to de- 
scribe the structure of *B, single particle potential models 
that treat ®B as a p-wave proton coupled to a ^Be core 
in its ground state, and microscopic cluster models that 
include two configurations of the three clusters ^He, a, 
and a proton. Recently, most experimenters have used 
the cluster model of Descouvemont and Baye (DB) [23 to 
extrapolate their data to zero energy. This generator co- 
ordinate method employs a central nucleon-nucleon inter- 
action along with Coulomb and spin-orbit interactions. 
Ref. [231 used the Volkov II nucleon-nucleon interaction 
|29j , while a more recent preliminary effort by Descouve- 
mont and Dufour (DD) ^] opts for the Minnesota force 
|3ll| . which describes low mass systems better. These 



cluster models include excited ^Be configurations. The 
single particle potential models generally employ a cen- 
tral Woods-Saxon -I- Coulomb potential, and optionally 
a spin-orbit interaction, which can be neglected in calcu- 
lations of the El S factor provided the central potential 
depth is properly adjusted. In these models, the p-wave 
potential depth is fixed by the *B binding energy. The 
depths for the other partial waves can be chosen identi- 
cally, but a better choice is to fix them using the well- 
measured s-wave scattering lengths for channel spin 1 
and 2 in the isospin mirror system ^Li + n |32| . The 
scattering lengths in the ^Be -I- p system have also been 
measured |33| , but with much lower precision. The scat- 
tering lengths for the dominant S ~ 2 channel are con- 
sistent between the isospin mirrors, but there is a 2.7a 
discrepancy in the value for the S = \ channel. This 
discrepancy is not understood at present, and deserves 
attention in the future. In this work, we consider the 
DD cluster model and the ''Li -I- n potential model of 
Davids and Typel (DT) [ll| , which reproduces the ^Li + 
n scattering lengths. These models well represent their 
respective classes. 

Both of these models have virtues. The cluster 
model allows more configurations than does the potential 
model, and therefore might be expected to describe the 
physics better. On the other hand, the potential model 
is simple, and has been tuned to reproduce the elastic 
scattering data. Fig. ^ shows the shapes predicted by 
the potential and cluster models. 



V. STRUCTURE-MODEL-DEPENDENT 
ANALYSIS 

With the theoretical structure models described briefly 
in Section llVl our fitting procedure involves only a single 
parameter, an absolute normalization. We fit data us- 
ing the DT ^Li -f n potential model and the DD cluster 
model employing the formalism presented in Section ITTl 
Our results are summarized in Tables III and Ull 

We first derive best fits for individual experiments. 
The 5*17(0) determinations are in excellent agreement 
with previous analyses [^[ij. In general, the discrepancy 
errors are smaller than or comparable to the intrinsic nor- 
malization errors of each data set. Note that the central 
values of 6*17(0) for the radiative capture data range from 
18-20 eV b and 21-22 eV b for non-Junghans and Jung- 
hans data sets respectively, while the CD data lie in the 
range of 16-19 eV b. This hints at some level of disagree- 
ment among the radiative capture data sets, especi ally 
when comparing the two Junghans experiments [8|, |20j 
with the other radiative capture data, in addition to that 
between the radiative capture and the CD data. 

To further explore and quantify this disagreement we 
look at different combinations of the data in our multiple 
experiment fits. We first look at the radiative capture 
data alone. In these model-dependent analyses, we con- 
sider 3 combinations of the radiative capture data, (1) 



TABLE I: Best fits to radiative capture (RC) and Coulomb 
dissociation (CD) data for Ecm < 425 keV using the DT ^Li ± 
n potential model. Shown are the best fit astrophysical S fac- 
tors and their standard deviations (statistical) at Ecn — 0. 
Also shown are the two individual contributions to the to- 
tal systematic error, the intrinsic normalization error and the 
discrepancy error, given in percent. An additional potential 
model parameter uncertainty should be added to the other 
systematic errors, with tmodei = 0.01%, reflecting the small 
uncertainties in the '^ Li + n elastic scattering data. 





Data set (# of points) 


Si7(0) (eV b) 


^norin 


^disc 




Filippone (6) 


19.3 ±0.5 


11.9% 


5.6% 




Strieder (2) 


17.9 ±0.6 


8.3% 


2.7% 




Hammache-1 (3) 


19.4 ±0.6 


4.9% 


5.7% 


Direct 


Hammache-2 (2) 


18.8 ±1.2 


12.2% 


8.1% 




Junghans-BEl (8) 


21.6 ±0.2 


2.7% 


0.7% 




Junghans-BE3 (13) 


21.2 ±0.1 


2.3% 


1.3% 




Baby (3) 


19.8 ±0.7 


2.2% 


2.2% 




Davids (2) 


16.6 ±0.5 


7.1% 


0.1% 


CD 


Schiimann (2) 


18.4 ±0.8 


5.6% 


0.8% 




All RC but Junghans 


19.3 ±0.7 


— 


3.8% 


indep. 


Junghans 


21.4 ±0.4 


— 


0.8% 


norm. 


All radiative capture 


20.8 ±0.4 


— 


5.1% 




Coulomb dissociation 


17.5 ±0.9 


— 


5.5% 




All but Junghans 


19.1 ±0.3 


8.5% 


6.2% 


corr. 


Junghans 


21.3 ±0.1 


2.6% 


1.4% 


norm. 


All radiative capture 


21.2 ±0.1 


5.6% 


3.2% 




Coulomb dissociation 


17.1 ±0.5 


6.7% 


5.5% 



TABLE H: Same as tableQ] except for the DD cluster model. 
A model uncertainty should be added to the other systematic 
errors, but DD do not present formal errors for this model. 





Data set 


5i7(0) (cV b) 


^norvn. 


^disc 


Direct 


Filippone 

Strieder 

Hammache-1 

Hammache-2 

Junghans-BEl 

Junghans-BE3 

Baby 


20.1 ±0.5 

18.8 ±0.7 
20.4 ±0.7 
19.1 ±1.2 
22.6 ±0.2 
22.1 ±0.1 

20.9 ±0.7 


11.9% 
8.3% 
4.9% 
12.2% 
2.7% 
2.3% 
2.2% 


6.1% 
2.4% 
5.9% 
8.2% 
1.0% 
1.4% 
2.7% 


CD 


Davids 
Schiimann 


17.4 ±0.6 
19.2 ±0.9 


7.1% 
5.6% 


1.0% 
1.6% 


indep. 
norm. 


All RC but Junghans 

Junghans 

All radiative capture 

Coulomb dissociation 


20.3 ±0.7 
22.3 ±0.5 
21.8 ±0.4 
18.2 ±1.0 


— 


4.2% 
1.1% 
4.9% 
5.3% 


corr. 
norm. 


All RC but Junghans 

Junghans 

All radiative capture 

Coulomb dissociation 


20.0 ±0.3 
22.2 ±0.1 

22.1 ±0.1 
17.9 ±0.5 


8.6% 

2.5% 
5.5% 
6.7% 


6.5% 
1.6% 
3.3% 

5.5% 
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FIG. 1: (Color online) Shapes of the S factors predicted by 
the Descouvemont and Baye (DB) and Descouvemont and 
Dufour (DD) cluster models and the Davids and Typel (DT) 
potential model as a function of relative energy (Ecm). Shown 
is the ratio of the El S factor to its value at Ecm = 0, which 
allows the shapes of models with different absolute normal- 
izations to be shown in the same figure. The shapes of the 
two cluster models are very similar at low energies, but sub- 
stantial deviations exist above 700 keV. The potential model 
is seen to have a slightly larger slope at high energies, but 
also a very different shape below 400 keV. 



all but the Junghans data, (2) only Junghans data, and 
(3) all the radiative capture data. Using the DT ^Li ± n 
potential model, the independent normalization method 
gives for 5i7(0) in cV b: (1) 19. 3 ±1.0, (2) 21.4 ±0.5 and 
(3) 20.8±1.1; using the correlated normalization method, 
we find (1) 19.1 ± 2.0, (2) 21.3 ± 0.7 and (3) 21.2 ± 1.4 
eV b. For the CD data we find 17.5 ± 1.3 and 17.1 ± 1.6 
eV b using the independent normalization and correlated 
normalization methods respectively. 

We find that the CD data and the non- Junghans radia- 
tive capture data are consistent, showing approximately 
a 1(7 discrepancy, while the CD data and Junghans data 
are discrepant at a little over the 2(T level. Furthermore 
there is a 1 to 2a disagreement between the Junghans 
data and the other radiative capture measurements. The 
Junghans data dominate the fit when combining all the 
radiative capture data, due to their extremely small er- 
rors and the size of the data sets. Similar but somewhat 
higher results are obtained with the DD cluster model. 
Interestingly, the high precision Junghans data are de- 
scribed slightly better by the DT ^Li + n potential model 
than by the DD cluster model, as seen in our quality-of-fit 
measure, the discrepancy error, but not at a significant 
level. The nature of the discrepancies between the Jung- 
hans data and both the non- Junghans radiative capture 
data and the CD data must be understood before we can 
address which model describes the data best. In an at- 
tempt to go beyond structure-model-dependent results. 



we now explore a structure- model- independent analysis. 



VI. POLE MODEL ANALYSIS 

Performing a structurc-modcl-independent analysis 
will provide insight into the quality of existing data for 
the ^Be(p,7)^B reaction below 425 keV. We use the ex- 
pansion suggested by [13, adopting a three-parameter 
model: 



Sn{E) ^ SniO) + a 



E 



QiE + Q) 



-f3E, 



(9) 



where Q = 137.5 keV. In addition to terms constant and 
linear in energy, this functional form contains a simple 
pole term, a universal feature of radiative capture reac- 
tions that is independent of the details of nuclear struc- 
ture models. We have chosen a fit that is linear in its 
parameters, but one can easily calculate the non-linear 
parameters adopted in jlOJ: the pole term is given by 
a = — a/S'i7(0) and the slope term by c = (3/Sn{0). 
A data set must have at least 3 data points in order 
for us to employ this three-parameter fit. Our results 
are summarized in Table IIIII Of all the data, only the 
most recent Junghans-BE3 data set |^ provides a signif- 
icant structure-model-independent constraint on the low 
energy behavior of the astrophysical S factor, implying 
517(0) = 24.3 ± 4.0 eV b. With this data set provid- 
ing the only significant low-energy constraint on Sn, it 
dominates the fit when we combine the radiative capture 
data to find 23.6 ± 3.4 eV b for the Junghans data and 
23.4 ± 3.4 eV b for all the radiative capture data, using 
the independent normalization method. Using the cor- 
related normalization method we find 24.1 ± 3.4 eV b 
for Junghans alone, and 18.9 ± 3.2 eV b if we include 
all data sets. With the current data we can place robust 
constraints on the extrapolated value <S'i7(0) independent 
of any particular structure model. We find that the cen- 
tral value is calculated to lie between 19-24 eV b with a 
total uncertainty of ±3.5 eV b. 

We see that in order to place a tighter constraint on 
S'17, one must assume some low-energy behavior. We now 
explore a constrained pole model which exploits knowl- 
edge of the astrophysical S factor around its subthreshold 
pole, on which all theories agree. 



VII. CONSTRAINED POLE MODEL ANALYSIS 

In order to pin down the S factor with higher precision, 
we need to examine the available theories and determine 
what, if any, universal information is available. To do 
this we fit each theory, normalized such that S'i7(0) = 1, 
using the pole model form, which now reduces to a 2 
parameter fit, depending on the pole term a = —a/SniO) 
and the slope term c = f]/ 817(0). We fit to 4 theories, 
the two potential models of DT and the cluster-model 



calculations of DB and DD. Our results are summarized 
in Table HVI 

As one can see the pole term is robustly determined, 
with a e [44,46] keV. We adopt a = 45 keV as our 
canonical value, and allow this information to propagate 
through our data analysis to see how our constraints on 
the low-energy behavior of S'17 improve. We thus use the 
two-parameter fit: 



Sn{E) = SuiO) 



1-a 



E 



Q{E + Q) 



-PE, 



(10) 



with the parameter a being fixed at 45 keV. Our results 
are summarized in Table Ivl 

Using this information, individual experiments deter- 
mine ^17(0) much more precisely than in the uncon- 
strained pole model fit. Again the Junghans data pro- 
vide the strongest constraints, 21.4 ± 1.3 eV b 2^ and 
21.4 ± 0.8 eV b j3] respectively. However, the discrep- 
ancy between the two Junghans data sets and the other 
radiative capture experiments is quite apparent, with the 
fit to the combined Junghans data giving 21.4 ±0.7 eV b 
and the other radiative capture data yielding 16.3 ± 2.4 
eV b. This is a 2a discrepancy between these two low- 
energy extrapolations, using the independent normaliza- 
tion method. This tension is somewhat reduced in the 
correlated normalization method, which yields 21.1 ± 0.8 
eV b and 17.9 ± 2.3 eV b for the Junghans and non- 
Junghans radiative capture data sets respectively, a dis- 
crepancy slightly more than Icr. The CD data yield 
17.5±2.5 eV b and 17.4±2.6 eV b using the independent 
normalization and correlated normalization methods, re- 
spectively. It is quite remarkable that the non-Junghans 
radiative capture data and the CD data agree so well, 
while both disagree with the Junghans data. This may 
be due to the sparseness of the data in these data sets, but 
the deviations from the Junghans results are significant. 
It would be desirable to have an independent confirma- 
tion of the Junghans measurements, since they deviate 
significantly from the other data and dominate the cen- 
tral value of combined fits. 

At this time, it is unclear what is causing the discrep- 
ancies. Arguably, the Coulomb dissociation measure- 
ments have the least in common with the radiative cap- 
ture measurements and large systematics of their own, 
so one could claim that they are the source of the dis- 
crepancy. However, this fails to explain the discrepancy 
among radiative capture measurements, namely between 
the two Junghans experiments [a, 123 and that of Filip- 
pone 15], which dominates the non-Junghans data sets. 
As shown in Table IVl none of the other radiative capture 
measurements below Ecm — 425 keV provides a signif- 
icant constraint on 6*17(0). It is unclear which experi- 
ments are responsible for the discrepancy. 

With these remaining uncertainties in mind, it is help- 
ful to remind ourselves that we have defined a rigorous 
treatment for exactly these kinds of discrepancies. Our 
treatment has examined the level of concordance and 
quantified it in terms of a discrepancy error. We thus 



recommend an astrophysical S factor at zero energy of: 

S'i7(0) = 20.7 ±1.2 eV b indep. normalization (11) 
5*17(0) = 21.4 ± 1.4 eV b corr. normalization (12) 

for the radiative capture data, and 

'S'i7(0) = 17.5 ± 2.5 eV b indep. normalization (13) 
817(0) = 17.4 ± 2.6 eV b corr. normalization (14) 

for the CD data. 

One can compare our low-energy 5* factor determina- 
tions with those determined from asymptotic normaliza- 
tion coefficients (ANCs). There are measurements from 
(1) proton transfer reactions [sj, (2) *B bre akup reac- 
tions [33 and (3) neutron transfer reactions |3a|. The 
asymptotic normalization coefficient can be very sim- 
ply related to the astrophysical 5" factor, so we quote 
the measurements in these terms. These determinations 
yield (1) 17.3 ± 1.8 eV b, (2) 18.7 ± 1.9 eV b and (3) 
17.6 ± 1.7 eV b for the astrophysical S factor of the 
^Be(p,7)®B reaction. These values of 5*17(0) agree per- 
fectly with the CD data and non-Junghans radiative cap- 
ture data. The ANC-derived 5* factors disagree with the 
Junghans data at slightly more than the la level. Again, 
these discrepancies need to be more fully explored to bet- 
ter determine the low-energy behavior of the S factor. 

Thus far we have primarily discussed the low-energy 
extrapolations of the 5 factor 5*17(0) and not the slope. 
Using our constrained pole model we find that the slopes 
based on the non-Junghans radiative capture, Junghans, 
and CD data are all roughly consistent with each other. 
In fact, the Junghans and CD slopes agree remarkable 
well, though the CD data have sizable errors. The non- 
Junghans radiative capture slope disagrees at the la 
level with both the Junghans and the CD data. Again, 
the non-Junghans fit is dominated by the Filippone [T^l 
data, with substantial uncertainties. These discrepan- 
cies disappear when one uses the correlated normaliza- 
tion method, suggesting that no significant deviation in 
the slope is observed and that the constraint on the slope 
using Ecm < 425 keV data is not particularly strong. The 
fact that the statistical errors in this parameter dominate 
over the normalization error supports this. 



VIII. CONCLUSIONS 

We have presented a robust formalism for fitting data 
that both properly propagates known systematic uncer- 
tainties and quantifies the quality of fit, incorporating 
a discrepancy error into the total systematic error. We 
discuss two limiting cases, one in which data sets are 
considered completely independent from each other, and 
a second in which all data sets are totally correlated. 
These two methods, the independent normalization and 
correlated normalization methods, yield similar results, 
providing robust constraints and suggesting that most 



previous analyses have underestimated the true uncer- 
tainty. 

A structure-model-dependent analysis was performed 
using the DT ^Li + n potential "ll] and the DD Min- 
nesota force cluster |3Qjl models. The ^Li -I- n poten- 
tial model generally predicts low-energy extrapolated 5* 
factors smaller than the Minnesota-interaction cluster 
model. With the available data, no significant prefer- 
ence for one model is observed. We explored a structure- 
model-independent fit to the data, finding that only the 
Junghans [g data placed even modest constraints on the 
low-energy 5* factor. Identifying a feature common to 
all models, the relative strength of the subthreshold pole 
term, we reduced the extrapolation error considerably. 
Even though we find evidence for discrepancies between 
the Junghans et al. radiative capture measurements and 
the others [isl \Wl \VA Ii§. |2]| ) , our rigorous and careful 
treatment of systematic errors provides a robust determi- 
nation of 5i7(0). Our analysis of indirect Coulomb disso- 
ciation data [nllli and ANC determinations |3ll3ll3^ 
found mutual agreement and consistency with the radia- 
tive capture measurements other than those of Junghans 
et al. 

The dominant source of error in the standard solar 
model predictions for the total ®B neutrino flux is the 
uncertainty in the heavy metal abundance. The domi- 
nant nuclear uncertainties stem from uncertainties in the 
■^He(a,7)^Be and ^Be(p,7)*B reactions. Using a tech- 
nique similar to that employed here, the authors of [3 
find a total error in the 5*34 normalization of 17%. If we 
adopt our determination of 5*17(0) using the independent 
normalization method (Eq. Ill|l and the 5*34 error assign- 
ment from 9] we find the following standard solar model 
(BP04) 6] prediction for the total ^B solar neutrino flux 



0(*B) = 5.63 [1 ± 0.058(5i7) ± 0.15(534) ± 0.21] . (15) 

Here we have separated the individual contributions to 
the total error in the neutrino flux, those from 5i7, 5*34, 
and the other standard solar model parameters, which 
when added in quadrature yield a total error of 26%. We 
can see that the new 5*17 error assignment contributes 
6% to the total neutrino flux error. 

While S17 now makes only a relatively small contri- 
bution to the total uncertainty in the predicted ^B solar 
neutrino flux, ongoing measurements of S34 and improved 
radiative opacity tables may reduce the other solar model 
uncertainties substantially in the near future. In order to 
further reduce the uncertainty on S17, a new high preci- 
sion measurement with a detailed error budget would be 
required. 
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TABLE III: Best fit parameters for data witfi Ecm < 425 keV witli tfie pole model functional form of Eq. |^ 





Data 


5i7(0) (eV b) 


a (eV b MeV) 


/3 (eV b MeV"^) 


^norm 


Iodise 




Fillipone 


38.6 ±15.7 


-8.5 ±5.4 


70.6 ± 37.2 


11.9% 


3.1% 




Hammache-l 


-2430 ± 2160 


-623 ± 554 


-2308 ± 2103 


4.9% 


0.0% 


RC 


Junghans-BEl 


18.4 ± 10.6 


-0.0 ±3.2 


7.1 ±17.5 


2.7% 


0.7% 




Junglians-BE3 


24.3 ±3.9 


-2.0 ±1.3 


18 ±9 


2.3% 


1.3% 




Baby 


55.3 ±213.1 


-11.4 ±56.2 


66 ± 229 


2.2% 


0.0% 




All RC but Junghans 


34.9 ±13.2 


-7.2 ±4.3 


60.2 ± 27.4 


— 


1.3% 


indep. 


Junghans 


23.6 ±3.4 


-1.7 ±1.1 


16.1 ±7.1 


— 


0.7% 


norm. 


All radiative capture 


23.4 ±3.3 


-1.8±1.1 


16.6 ±6.8 


— 


3.7% 




Coulomb dissociation 


— 


— 


— 


— 


— 




All RC but Junghans " 


35.4 ± 13.0 


-7.1 ±4.2 


57.2 ±26.8 


6.7% 


3.6% 




All but Junghans ' 


29.0 ±7.9 


-5.0 ±2.9 


43.5 ±21.2 


7.7% 


5.4% 


corr. 


Junghans 


24.1 ±3.4 


-1.9±1.1 


18.5 ±6.9 


2.6% 


1.3% 


norm. 


All radiative capture 


18.9 ±2.9 


-0.1 ±1.0 


5.7 ±6.1 


6.0% 


3.2% 




Coulomb dissociation 


-13.5 ±26.4 


9.3 ±8.6 


-53 ± 53 


6.5% 


3.6% 



"All non-Junghans RC data sets with at least 3 points 
''All non-Junghans RC 



TABLE IV: Best fit parameters for several theoretical struc- 
ture models below Ecm ~ 425 keV. 



Model 


a (keV) 


c (MeV"^) 


DT ^Li ± n potential 


45.7 ±0.5 


0.553 ±0.007 


DT ''Be ± p potential 


45.0 ± 0.6 


0.433 ± 0.007 


DB Volkov II 


45.5 ±0.1 


0.434 ± 0.001 


DD Minnesota 


44.5 ±0.1 


0.404 ± 0.002 
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TABLE V: Results of a constrained pole model fit of radiative capture and Coulomb dissociation data below Ecm ~ 425 keV. 
Shown are the best fit astrophysical S factors, the associated slope parameter /3, and their standard deviations (statistical). 
Also shown are both the intrinsic normalization errors and the discrepancy normalization errors for combined data sets, all 
cited in percent, which should be added in quadrature with the statistical errors to get the total error. An additional systematic 
error must be added due to the uncertainty in a; ASn{0)/Sn{0) ~ 0.16Aa/a, and A/3//3 « 0.80Aa/a, which for a = 45 ± 1 
keV yields A5'i7(0)/5'i7(0) = 0.4% and A/3//3 = 1.8%, respectively. 





Data set 


5i7(0) (cV b) 


P (eV b MeV"^) 


^norm 


^disc 




Filippone 


16.4 ±2.8 


18.1 ±7.2 


11.9% 


4.9% 




Strieder 


30.6 ± 17.2 


-17.8 ±37.6 


8.3% 


4.0% 




Hammache-1 


-2.7 ±16.7 


56.2 ±34.4 


4.9% 


3.7% 


RC 


Hammache-2 


-2.3 ±16.1 


159. ±113. 


12.2% 


0.0% 




Junghans-BEl 


21.4 ±1.1 


12.1 ±2.6 


2.7% 


0.7% 




Junghans-BE3 


21.4 ±0.6 


11.3 ± 1.6 


2.3% 


1.3% 




Baby 


14.6 ±8.5 


21.9 ±17.8 


2.2% 


0.7% 




Davids 


16.6 ±2.5 


9.1 ±6.4 


7.1% 


0.0% 


CD 


Schumann 


17.5 ±4.9 


12.5 ±13.2 


5.6% 


0.0% 




All RC but Junghans 


16.3 ±2.3 


17.3 ±5.0 


— 


3.5% 


indep. 


Junghans 


21.4 ±0.7 


11.6 ± 1.4 


— 


0.8% 


norm. 


All radiative capture 


20.7 ±0.6 


11.3 ±1.3 


— 


4.8% 




Coulomb dissociation 


17.5 ±2.3 


9.7 ±5.8 


— 


5.1% 




All RC but Junghans 


17.9 ±1.4 


13.4 ±3.3 


8.0% 


5.8% 


corr. 


Junghans 


21.1 ±0.5 


12.3 ±1.3 


2.6% 


1.3% 


norm. 


All radiative capture 


21.4 ±0.5 


11.0 ± 1.2 


5.7% 


3.2% 




Coulomb dissociation 


17.4 ±2.2 


8.7 ±5.7 


6.4% 


4.8% 



